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INTRODUCTION 


• IcIk??,o£2 SITI0N J0,NT PDF METHOD HAS 
BEEN USED TO MODEL A WIDE CLASS OF 

GASEOUS TURBULENT REACTIVE FLOWS. 
(S.B. POPE) 

• NONLINEAR CHEMICAL REACTION RATES 
COULD BE EVALUATED WITHOUT ANY 
APPROXIMATION. 

• AN EXTENSION OF THE PDF METHOD TO THE 
MODELING OF SPRAY FLAMES. 

• EVALUATE THE LIMITATIONS AND 
CAPABILITIES OF THIS METHOD -IN THE 
MODELING OF GAS- TURBINE COMBUSTOR 


Composition Joint Pdf Transport Equation 


PP.t + puip^ + [pw a (tp)p]^ = 

{Mean convection } {Chemical reactions} 

~\P < < I ± > pU - [p < I ± > pU„ 

{Turbulent convection} {Molecular mixing} 
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Density- weighted joint pdf. 
chemical source term for the a-th 
composition variable, 
conditional average of Favre velocity 
fluctuations. 

condi tional average of scalar dissipation, 
conditional average of liquid-phase source 
term for the a-th composition variable. 
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Modeling Aspects of the Pdf Transport Equation 


• < u" | ?/> > is modeled using a gradient-diffusion model. 

t < I | > is modeled using a variant of Curl s model. 

• The new term <’S a \}P> involving the conditional average 
of liquid-phase source term is modeled based on the average 
values of species and enthalpy: 

< \ ■>!>>= TT— ~ M 

p ~ P&v i -' 

for <j> a s= Y a % Of = U 2, ...,s = (T — 1 

< i.S a | >= n k m k(~ l k.cfJ + - &*) 
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MODELING ASPECTS 


THE MODELED PDF TRANSPORT EQUATION 
PROVIDES THE SOLUTION FOR JH| SPECIES 
AND TEMPERATURE FIELDS WITH THE MEAN 
VELOCITY AND THE TURBULENT 
DIFFUSIVITY AND FREQUENCY PROVIDED AS 
INPUTS FROM THE CFD SOLVER AND THE 
SPRAY SOURCE TERMS FROM THE 
LIQUID-PHASE SOLVER. 

THE MEAN FLOW AND TURBULENCE 
EQUATIONS ARE SOLVED BY A 
CONVENTIONAL CFD SOLVER WITH THE 
MEAN SPECIES AND TEMPERATURE FIELDS 
PROVIDED AS INPUTS FROM THE PDF 
SOLVER AND THE SPRAY SOURCE TERMS 
FROM THE LIQUID-PHASE SOLVER. 

THE LIQUID-PHASE EQUATIONS ARE 
FORMULATED IN a ^GRANGIAN 

COORDINATES WITH APPROPRIATE 
CONSIDERATION TAKEN INTO ACCOUNT OF 
THE EXCHANGES OF MASS, MOMENTUM, 
AND ENERGY BETWEEN THE TWO PHASES. 
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NUMERICAL METHOD 


• Mean-Flow and Turbulence Equations 

- Axisymmetric, Unsteady. 

- Incompressible Navier-Stokes (Variable-Density). 

- A Standard Two-Equation k-e Turbulence Model. 

- A Pressure-Based CFD Solver Based on the 
SIMPLE Algorithm of Patankar and Spalding. 

• Liquid-phase Equations 

- The Spray Model (Raju and Sirignano). 

- Dilute Spray Assumption. 

- The ODE’s for the Particle Size, Velocity, and 
Location are Solved Using a R-K Method. 

- The PDE’s for the Internal Droplet Distribution 
(Vortex Model) are Solved by an Implicit Method. 

- Droplet Regression Rate is Based on Either a 
Gas-Phase Boundary Layer-Analysis or 
Low-Reynolds Correlation. 


NUMERICAL METHOD 


• The PDF Transport Equation 

- A Fractional Step Monte-Carlo Method (Pope). 

- Spatial Transport, Molecular Mixing, Liquid-Phase 
Source Terms, and Chemical Kinetics are 
advanced in a Series of Sequential Steps. 

- Vectorization 


• Interaction Between the Two Phases 

- Interpolation of the Gas-Phase Properties at the 
Particle Location Using an Area-Weighted 
Averaging. 

- The Source Terms Evaluated at the Particle 
Location are redistributed among the surrounding 
Computational Nodes Using an Area-Weighted 
Averaging. 
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CHEMICAL KINETICS MODEL 


• IT IS BASED ON A SINGLE STEP GLOBAL 
MECHANISM OF WESTBROOK AND DRYER 
FOR N-DECANE/OXYGEN COMBUSTION. 


• THIS GLOBAL COMBUSTION MECHANISM 
WAS SHOWN TO PROVIDE ADEQUATE 
REPRESENTATION OF TEMPERATURE 
HISTORIES IN FLOWS NOT DOMINATED BY 
LONG IGNITION DELAY TIMES. 



Geometry of the combustion chamber. 
(El Banhawy and Whitelaw) 
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EXPERIMENTAL DETAILS 


• The experimental data corresponds to 
the following inflow conditions: 

inflow temperature = 31 0 K, 
air mass flow rate = 355 kg/h, 
air/fuel ratio =20.17, 
swirl vane angle = 45 deg, 
swirl number =0.721. 

• The reported error in the 
measurements is about 10 to 15 % for 
the temperature and about 15% for the 
velocity. 


Details of Fuel Injection 

• A fuel nozzle of swirl-atomization type was used. 

• The liquid fuel injection is simulated by injecting 
a discretized parcel of liquid mass at the end of 

each A t injection 

• The droplet-size distribution is given by: 

£- 4.21 10 « 

n l A )jJ Am 



Drop 1 1 1 dldmater, nlcrons 


Droplet size Distribution 

• The initial droplet injection velocity corresponds 
to: Uk — 11.0 m/s, Wk — 6.1, and Vk — 0.5 — 2.5. 
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PARAMETER SELECTIONS 


• The computations were performed on 
a grid with a mesh size of 60x60. 

• The PDF solution is obtained by 
making use of 250 particles per cell. 

• Dtg - Dtj n j ec tion = 1-5 ms, Dt^ = 
0.0375 ms, and Dt Monte _ Car|o = 
0.015 ms. 

• Two CPU seconds on a CRAY Y-MP 
per one Dtg and about 2 to 3 CPU 

hours - 4000 time steps - for the 
solution to reach steady state. 



Temperature contours and droplet locations. 
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Near wake radial profiles of velocity. 



Schematic of an open spray burner. 
(Dan Bulzan of IFMD at LeRC) 
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CONCLUDING REMARKS 


• The comparisons show that the general features of 
the flowfield are correctly predicted by the present 
solution procedure. 

• The present solution appears to provide a better 
representation of the temperature field, particularly, in 
the reverse-velocity zone. 

• The overpredictions in the centerline velocity could be 
attributed to the following reasons: 

- The use of k-e turbulence model is known to be less 
precise in highly swirling flows. 

- The swirl number used here is reported to be 
estimated rather than measured. 
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